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Electron counting experiments attempt to provide a current of a known number of electrons per unit time. We 
propose architectures utilizing a few readily available electron-pumps or turnstiles with modest error rates of 1 
part per 10 4 with common sensitive electrometers to achieve the desirable accuracy of 1 part in 10 8 . This is 
achieved not by counting all transferred electrons but by counting only the errors of individual devices; these are 
less frequent and therefore readily recognized and accounted for. Our proposal thereby eases the route towards 
quantum based standards for current and capacitance. 



I. INTRODUCTION 

The long pursued goal of closing the quantum metrological 
triangle HKH involves the comparison of a current given by 
Josephson and Quantum-Hall effect, with one based on count- 
ing electrons flowing through a constraint. An equivalent al- 
ternative is to charge a capacitor with a counted number of 
electrons J4). 

For metrologically precise measurements of current about 
1 nA is required [5|. This can be provided either by pump- 
ing at frequencies higher than those that have shown satis- 
factory accuracy or by parallelization. Higher operation fre- 
quencies increase the error-rates, while parallelization is hin- 
dered by conflicting requirements on the devices, including 
yield and accuracy. We propose here a concept to alleviate 
both concerns; this is made possible by the relative simplicity 
of non-adiabatic single gate devices providing quantized cur- 
rent, which have been realized in several technologies |6H9) 
and even parallelized iflOl [TTI 

On the other hand error-rates akin to those of the shut- 
tling experiment conducted for adiabatic pumps |4] |3 have 
not been demonstrated, while direct current measurements in- 
dicate that error rates near T = 1 ppm can be achieved ifTZl . 

To improve upon these results, we propose a circuit that 
utilizes the high precision and accuracy provided by electron 
pumps, but permits the correction of the remaining errors: We 
connect several pumps in series and observe the charges on 
memory nodes between them. An example of such a circuit is 
shown in fig (1). 

Each of the aforementioned pumps will provide a current 
very close to eQ, when its gate is driven with periodicity I/O, 
— e being the electron charge; and will barely pass electrons 
between the terminals in a suitable off state. 
When all turnstiles move electrons with the same rate, then 
only each node-charge averaged over a pump-period is con- 
stant. Whenever an error, a surplus (or deficit) electron being 
shifted by pump j, occurs, then the charge on node j — 1 will 
drop by one electron, while the charge on node j will rise by 
it, labeling the outside reservoirs as nodes and N. Accord- 
ingly, errors committed by different pumps result in different 
charge signals on the nodes. Such different signatures can be 
recognized and attributed. If this were to be done flawlessly, 
it would be possible to tell with certainty by how many elec- 



trons the charge transported by each pump deviates from the 
intended value. 

But when the characteristic time to identify the charge state, 
tm, and the time between errors, 1/Tfi, are comparable then, 
if all observed errors are attributed to the pump (or pumps) 
whose error would most likely explain the observed signature, 
some multi-pump errors will mistakenly be mis -attributed: 
Events, where the first two pumps of a three pump circuit 
both pump one electron more than intended, cannot be distin- 
guished from the -more likely- case of one too few electrons 
pumped by the third pump. Only a four stage circuit would be 
able to distinguish these two scenarios. Circuits with an odd 
number of pumps AT exhibit 2N\/((N - 1)/2)!((A + l)/2)! 
different scenarios of errors of (N + l)/2 pumps that can- 
not be distinguished from errors of all other pumps. If these 
happen within a time-interval tm, then they will be mis- 
attributed. The relative rate of such mis-attributions is then 
approximately 

^corrected — 2 N _ 1 | N+1 | (rOTjVf ) 2 Qrj^ ' ^ 

So the corrected error-rate scales as the bare error rate, T, 
risen to the power of half the number of stages employed; this 
allows devices with a sufficiently low bare error-rate to reach 
any desired corrected error-rate with additional stages; this 
assertion is demonstrated numerically in section (III) of this 
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FIG. 1 : Schematic of a circuit for an error-correction scheme in a typ- 
ical electron counting experiment composed of metallic single elec- 
tron transistors and electron pumps utilizing dynamic quantum dots 
in a 2-dimensional electron gas. 
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paper after the required analytical framework is provided in 
section (II). 



II. MASTER EQUATION TREATMENT OF THE 
ARBITRARY ERROR-CORRECTION CIRCUIT 

We identify the net errors the N pumps have committed in 
time interval [0, t) with the stochastic process & with outcome 
space Z N over time and denote the accompanying probability 
distribution over Z N with <p(t). With i = i 2 , ...i N ) the 
entries of vector <f> are labeled <fi(i)(t) and are the joint prob- 
abilities of ij errors having been being committed by pump 
j for any j in the time-interval [0, t). As the error-events we 
intend to observe are assumed to be rare, we can assume K 
to be a Markov-process. For the discrete time series at inte- 
ger multiples of a small time St process 8. is described by the 
Pauli-Master equation 

4>{t + St) - ${t) = K$(t)St, (2) 

where K is the stochastic matrix describing J? and constant 
if 8. is stationary. K is composed of the probabilities of an 
error F occurring during a time-interval [kSt, (k + l)St) given 
initial state of the errors i as Pp(kSt,i) — StK^ + p^. The 
probability for the realization of a particular trajectory F, with 
F(kSt) errors committed in each time interval [kSt, (k + l)St) 
is then 

T/St fe-1 
fe=0 k'=0 

In addition to the random-walk of eq. (|3]l we utilize a number 
s of noisy measurements in the circuits of interest. For each 
time the outcome of the s measurements can be written as a 
vector £ in W. We assume that, given state i of the random- 
walk the measurement outcomes follow the form 

(kSt) = f{i(kSt) , kSt) + n(kSt) , (4) 

where n(t) is stochastic additive noise and / the determinis- 
tic transfer-function relating observed system to measurement 
outcomes. We assume the noise to be independent of and un- 
correlated to the random walk, to be stationary, ergodic and 
markov and, crucially, to be uncorrelated in time, which is 
a valid approximation, when the errors are rare compared to 
the noise's auto-correlation time. For the noise averaged over 
interval of length St we denote the noise probability-density 
over W with vgt(t, n). Equivalently to the time-trajectories of 
the error-scenarios T, we define S to be time-trajectories of 
the measurement outcomes and note that the probability that 
the circuit commits a particular random walk T and generates 
electrometer outputs S over interval [0, T) is 

T/St k 

nf = II MkSt,S(kSt) - f{J2 F{k5t),kSt)), (5) 

fe=0 fe'=0 



By Bayes theorem the conditional probability IT^s of a given 
random walk T having occurred when output S is exhibited 
is 

n^ 5 = nf/^nf, (6) 

whenever the denominator is non-zero. Accordingly, given S 
(over the entire time interval [0, T)) the probability that net 
errors E have occurred in the time-interval [0, t) is the sum of 
these conditional probabilities over all possible paths T satis- 
fying said condition : 

n(E,t\s)= nf/£nf . (?) 

While the probability of E errors having occurred in time in- 
terval [0, t), given measurement S over that time interval and 
any measurement outside of it, is 

n(g,t|5,o = E y' ) ^ (t,)vt ' £tS ; l ^" r . (8) 

l^S'(t')=S(t')\/t'<t L-<T 

Where the interest is not in identifying one particular random 
walk per se, but to identify the probabilities that the final result 
of the underlying random walk has been E net-errors, pro- 
vided measurements S; in particular if there is no need (or 
possibility) to identify in what precise time interval an error 
occurred, then indeed H(E,T\S ,T) is going to be the vari- 
able of interest. Obviously no more can be desired, once St is 
chosen sufficiently small. 

More convenient to compute is the probability that up to time t 
E net errors have occurred and measurement-outcome S been 
observed over the same interval, but irrespective of S for times 

t' > t: 

m,t\s,t)= Yl E n 5 ; - w 

S'(t')=S(t')Vt'<t F^E.t 

When stepping forward in time by St, the change of 
ip(E, t\S, t) then includes four terms: random walks leaving 
the second sum, because of an error-event, F{t + St) ^ 0; 
random walks entering the second sum for the same reason, 
F(t + St) ^ and J2t'<t+St ^(*) = ^ and ^'(0 = 
S(t')yt' < t + dt; and paths S' leaving the first sum, be- 
cause S'(t + St) ^ S(t + St), but satisfying F(t + St) = 0, 
to avoid double counting with respect to the first term. Of 
course there are as well paths that stay within both sums, with 

E*<t = Et/< t +«B p (t) = £ and s'V) = £(*')vf < 

t + St; while no paths can enter the first sum by its definition. 
Accordingly 
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tp(E, (k + l)St\S, (k + l)St) - tp(E, kSt\S, kSt) = - ^ Pp(kSt, E)ip(E, kSt\S, kSt) 

F{kSt)^0 

+v(S((k + l)St) - f(E, (k + l)St)) Pjp(k5t,E-F)ip(E-F,k5t\S,kSt) 

F(kSt)jtO 

-[1 - u{S((k + l)5t) - f(E, (k + l)5t))]P^(E,t\S,t) 



(10) 



where is the null-vector. For each time ip contains for 
each trajectory S a number of elements equal to the num- 
ber of possible net error scenarios E. For a given trajectory 
with N pumps under consideration, at each time the set of all 
ip(E, t, S, t) will be an element of [0, 1] Analogous to </> it is 
convenient to summarize all ip(E, t, S, i) into i/j(t). Dropping 
the reference to S. Eq. ( [T0| simplifies to 



tory S is then computed as 

<E\s>=y: 

* ms = Y^{E-<E\S>Y *&T> S ,T) 



4>{t + 8t) = (l + USt)(l + Kdt)ip(t) 



e u5t (l 



■K6t)ip{t). 
(11) 

Here U5t includes the dependence on measurements v of eq. 



( 10 1 and is a diagonal matrix; while A'does not depend on S 



and is constant if the error-rates do not depend on time, that is 
if & is stationary. Remarkably the evolution of tp is governed 
by a linear difference-equation akin to eq. pi. While the defi- 
nitions both Pp and v depend on the choice of St, we assume 
that it can be chosen sufficiently small, such that neither K 
nor U do. 

For all times t > 0, it is obvious that 



E S'(i')=S(t')W<t F 

(12) 

so it is convenient to define that 



^(0, 0,5,0) = ££V(£,0,S,0) = 1, (13) 



Given measurement outcome S this computes that (on 
average) < E\S > net-errors have occurred in time-interval 
[0, T]; that quantity has an uncertainty of (Jg^ s , however the 

entire probability distribution tp(T) is known. The mean of 
that latter quantity weighted over all possible trajectories will 
be of particular interest below, < crgi s >s~ J2s^ s<7 e s >^ e 
expectation-value of the uncertainty of the error-accounting 
scheme. For computer-generated S and T the average error 
of the error-accounting scheme can be computed as well, as 
the difference, 5%, of the computed < E\S > and ex-ante 
known underlying error of the random walk. However, 
< fxg| 5 >s is a sufficient proxy, that is accessible when only 
S is known, as in an experiment. 

For any random walk governed by a known stochastic matrix 
K and additive totally-random noise described by U, eq. 



(Ill is readily integrated for any one set of measurement 
outcomes. Expectation- value of the occurred error and its 
uncertainty according to eq. ( fT~4] > follow. Average over all 
possible outcomes, however, can only be approximated by 
summing over a large set of trajectories. 



which implies that the set of all possible measurement out- 
comes S prior to the first measurement has one element of 
probability 1. With initial condition eq. ( 13 1 and time devel- 



opment of eq. ( 10 1 or eq 
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i the final probability of E errors 
having occurred in interval [0, T) and trajectories S having 
been measured can readily be computed as if>(E, T, S, T); the 
conditional probability of E errors having occurred in inter- 
val [0, T) given trajectories S having been measured is then 
tp(E, T, S, T) I J2p' ^(E 1 , T, S, T), where the denominator 
is just the 1-norm of vector 4>(T)- Eq. (Ill is the master- 



equation governing the net-error probabilities of the shrinking 
set satisfying trajectories S. 

The expectation value and variance of the errors given trajec- 



III. NUMERICAL EVALUATION OF THE 
ONE-DIMENSIONAL PUMP CHAIN 

The analytical treatment above is valid for an arbitrary cir- 
cuit, while for numerical studies we confine ourselves to the 
simplest case, which is of particular experimental interest, the 
1 -dimensional chain of electron pumps. In eq. (Hh this corre- 



sponds to fj(i(t), t) 



ij+i, where subindex j indicates 



the j-th component of the vector of length s = N — 1. The 
electrometer noise we assume to be ergodic, stationary, gaus- 
sian, white. The random walk we assume to be stationary - 
though pulsed operation of the pumps may be of interest - and 
independent of the charges accumulated on the nodes. This 
case is understood even using only the arguments yielding eq. 
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FIG. 2: Signal traces of first (black) and second (gray) electrometer. 
Black and red lines indicate ideal measurement results without noise. 
Pump-error rates are 1000/s, per pump and sign, electrometer noise 
10 _5 e/\/Hz, coupling capacitance ratio Cc/Cnode = 0.04. 



FIG. 3: Corrected error rates with respect to time for three (upper 
curve) and five (lower curve) pump circuits as a function of the un- 
corrected error rates. Dashed lines are a fit to eq. {JJ. Other parame- 
ters as in fig. [2] 



([TJ and could experimentally be realized by applying feed- 
back voltages to the nodes. Then the net-errors will diffuse 
freely for each pump and will accordingly diverge with y/T. 
As the electrometers cannot distinguish states that differ by 
the same number of errors for all pumps, we combine charge 
equivalent states when computing ip. For free error diffusion 
we need to compute ijj only over a few states around the most 
likely error-scenario, typically all states that deviate by up to 
three errors from the most likely error-scenario. Integration 
of eq. ( [Tol l i s simplified as K is both constant and invari- 
ant under translation by any possible error. The latter is no 
longer the case when voltage dependence of the error-rates is 
included; the constant K is here only invariant under uniform 
translations of errors on all pumps. Accordingly, dt/j/dt needs 
to be calculated over a larger configuration space. This is the 
case of mesoscopic feedback fT3l [Mil ; however such charge 
dependence of the error-rates does not add conceptual compli- 
cation to our treatment. Figure [2] shows a typical set of signal 
traces that can be expected from the combination of the diffu- 
sion dynamics and electrometer noise realistic for the RF-SET 
ifTsHTT l and similar devices [ 18-20 1. These traces are used to 
numerically integrate eq. (Hi yielding a perceived diffusion. 



The result of this perceived diffusion is subtracted from the 
underlying diffusion. This difference is the error that our algo- 
rithm mis-attributes. Nonetheless, the rate of such errors after 
correction is much smaller than the rate of errors of each sin- 
gle device. Numerical results for rates of these reduced errors 
are shown for various parameters in fig. [3]in excellent agree- 
ment with eq. ([JJ). Even for a rate of 200 errors per second a 
1000-fold improvement of the error rate is achieved with only 
three pumps and two intermediate electrometers, so that even 
devices exhibiting a relative error-rate of 10 ppm can reach 
the metrologically desirable error threshold of 10 ppb, due to 
the proposed architecture. The five pump version achieves 
a 10000-fold improvement even at a rate of 3000 errors per 
second. Pumps with intrinsic error-rates of 100 ppm at pump 
rates of 30 MHz could thus be corrected to the required 10 ppb 
accuracy. If devices with 10 ppm were to be used, a 1000-fold 
improvement were to suffice, so that more than 8000 errors 
per second could be accounted for. This would correspond to 



a pump rate of 800 MHz, so that both the requirements of cur- 
rent magnitude and accuracy are met. Recall that the lowest 
reported error-rate for single parameter pumps is about 1 ppm 

m. 



CONCLUSION 

In summary, we have proposed a family of architectures 
utilizing either conventional electron pumps of modest qual- 
ity or one of the single-gate devices that are capable of pro- 
ducing quantized current steps with the now common RF- 
SET or a similar device. The algorithm constructed in this 
work determines the charge transported by such structure to a 
much greater precision than would be feasible by any single 
of the devices used. The relaxation on the requirements for 
the quality of the individual devices used, simplifies experi- 
mentation, even at higher operation frequency, and increases 
sample yield. It thereby opens the prospects towards paral- 
lelization needed to increase the provided current. We note 
that the proposed architecture allows to monitor the error-rates 
of individual pumps on the level of single electrons, similar 
to those stated for adiabatic electron-pumps in their shuttle 
mode [4|. This has not been possible for the single gate de- 
vices J6] [8] |9) in currently employed circuits. Accordingly 
the architecture proposed here allows further research on these 
devices and thereby a determination regarding their suitabil- 
ity for metrological applications, while easing the required 
thresholds. During the preparation of this manuscript we par- 
ticipated in work demonstrating the integrated operation of 
the components required for the experimental realization of 
this concept J2T|. 
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